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Abstract 

We show that the supersymmetric structure (in the sense of supersymmetric quantum 
mechanics) appears in Helmholtz optics describing light propagation in waveguides. For the 
case of elliptical waveguides, with the accuracy of paraxial approximation it admits a simple 
physical interpretation. The supersymmetry connects light beams of different colors. The 
difference in light frequencies for the supersymmetric beams is determined by the transverse 
gradient of the refractive index. These beams have the same wavelength in the propagation 
direction and can form a stable interference pattern. 


1 Introduction 

There is a correspondence between quantum mechanics and optics in the paraxial regime [1] 
and in the global one [2]. Consequently, many notions can be transported from one to the other, 
such as coherent and squeezed states (see, t.g n Ref. [3, 4]). Here we describe an optical system 
that exhibits supersymmetry in the sense of supersymmetric quantum mechanics (SUSY QM). 

Supersymmetry in quantum mephanics connects two Hamiltonians with the same spectrum 
except for the ground state (unbroken supersymmetry [5], see [6] and references therein.) SUSY 
has been successfully applied in atomic [7] and nuclear physics [8]. The supersymmetric form of 
the Dirac equation in an external field [9, 10, 11] presents the Dirac equation as the square root 
of the Klein- Gordon equation; the kind of supersymmetry we have here is analogous and finds 
the square root of the Helmholtz equation. The supersymmetric structure of Helmholtz optics 
describes light propagation in a waveguide and admits a very clear physical interpretation. 

We consider optical waveguides along the z-axis, i.e., media that are inhomogeneous only in 
the x direction. From the wave equation in 2 — F 1 dimensions for solutions of time frequency v we 
have the Helmholtz equation 

[<?* + d\ + i/ 2 n 2 (z)/c 2 ] f(x,z) - 0. (1) 

In Section 2 we present the supersymmetric structure of Eq. (1), in Section 3 a physical 
reinterpretation, and some concluding remarks in Section 4. Note that until Section 3 we do not 
use the paraxial approximation. 
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2 Supersymmetric structure 


We start with a system of two first-order equations for a two-component wave function of the 
x and z-coordinates, 


«(!;)-(: -S)(t)- 


(2) 


d 2 


^2 


' ) ( ) 

V.V+ ) \ ) 


where k is a constant and 

V± = ±d x + W(x), 

with an arbitrary function W{x). The second z-derivative involves the square of this 2x2 matrix, 

( —k 2 + v + v- 0 

V o - k 2 + 

Therefore we have two different Helmholtz equations for the components: 

[d 2 z + d 2 x + k 2 -W x - W 2 ] = 0, 

[d 2 z +d 2 x + k 2 + w x - W 2 ] = 0. (3) 

where W x = dWjdx. Now, for a given wavenumber k in z direction that is common to both 
wavefunctions, we write 

»i,a(*,*) = (4) 

Equations (3) then become the eigenvalue equations for the components, 


$1 = 

[/ fc 2 - K 2 } 

*1, 


]*a = 

[k 2 - K 2 

]*2. 

(5) 


Now we introduce the supersymmetric structure by considering the component <!>! as represent- 
ing the ‘boson’ and $2 the ‘fermion’ sectors of the supersymmetric Hamiltonian eigenfunctions. 
The supercharge operators are 


Q- = 


0 0 

t> + 0 


Q+ 


-c '•) 


The supersymmetric Hamiltonian, defined as the anticommutator of the two supercharges, is 


h, = = 


v_v + 0 
0 v+v 


.)• 


[H.,Q±] = o. 


The eigenvalue equations for H e are then 

w_v + 4>j = ( k 2 — /c 2 ) $ 1 , 

u + u _$ 2 = ( k 2 — /c 2 ) $ 2 - 

These equations coincide with Eqs. (5). Operators and v+v- have the same spectrum 

except for the ground state, which has to be a normalizable solution of one of the two equations 
v+<j>\ = 0 or v_<f >2 = 0. Therefore, SUSY connects solutions of Helmholtz equations of two different 
waveguides with index profiles 

v 2 n 2 (x)/c 2 = k 2 W x — W 2 . (6) 
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3 Physical Reinterpretation 


We now reinterpret these formulas in another way to describe two different beams in the srnne 
waveguide. We choose 

W{x) = wx. (7) 

The eigenvalue equations (5) are then two Schrodinger equations for harmonic oscillators with 
displaced energy levels, 

[-d 2 + w 2 x 2 ] $i |2 = (A: 2 =F uj - k 2 ) $i, 2 - ( 8 ) 

Meanwhile, seen as Helmholtz equations, Eqs. (8) correspond to 

u 2 n 2 (x)jc 2 = k 2 ^w- w 2 x 2 , 


( 9 ) 


with the same refractive index n(x). In the absence of material dispersion, n is independent of 
v. This is an elliptic index-profile waveguide; it approximates the usual parabolic index-profile 
waveguide in the paraxial regime [12]. The eigenvalues of the operator on the left of (8) are well 
known to be 

E m = w(2m + l), m = 0,1,2,.... (10) 

We now find conditions for a supersymmetric pair of Helmholtz solutions propagating in the 
same waveguide. Into this waveguide n 2 (x ) = tIq — n 2 x 2 we inject two light beams with slightly 
different time frequencies 

<2 = ^o(l T c )- 

To fulfill Eqs. (9), we substitute k — n 0 L> 0 /c } w — n\V$jc and find the conditions 

£ = ^ = ^ = 2£ - (H) 

This determines the frequency shift in terms of the transverse index profile of the waveguide. 

From the standard quantum harmonic oscillator solution we know that the Gaussian beam 
width in x-direction is 'Ax = (2tu) -1 ^ 2 . Therefore, the parameter 

e — [2ib 2 (Ax) 2 ] -1 ~ (wavelength/beam width) 2 

should be small. The two light beams with the frequencies v h2 obey the Helmholtz equations with 

z/ 2 2 n 2 (x)/c 2 = k 2 w - w 2 x 2 (l T e). (12) 

When c — > 0, these two Helmholtz equations become a supersymmetric pair that is slightly broken 

by the term tw 2 x 2 . * i . . 

It is easy to see that SUSY has the same accuracy as the paraxial approximation. Solution 
of the Helmholtz equation (1) with the term (12) (compare with Eqs. (8)) shows that the z- 
propagation wavenumbers of the two beams are 


/c 2 >2 (m) = k 2 - w[( 1 ± e)(2m + 1) ± 1], m = 0, 1, . . . 


or 


Ki >2 (m) 


~ 1 — e ( m 
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( m+ 2 ± 2. 

) 1 

l m + 2) 


+ 


377 



The term of order e — w/k 2 exhibits exact SUSY; except for the ground state k 2 (0) = k, the 
z-wavenumbers of the two beams coincide 

Ki(m - 1) = K 2 (m) + 0(e 2 ), m = 1,2,*-*, (13) 

The two terms of order e 2 respectively give a nonlinear correction to the paraxial approximation 
and break supersymmetry. Therefore, SUSY is exact in the paraxial approximation and is broken 
beyond this regime. 

Thus, supersymmetry (13) connects light beams of different frequencies v 2 2 in the same wave- 
guide [cf. Eq. (11)], but having the same wavelength 2tt/k in the propagation direction z [Eq. 

(4)]. These two beams form a stable interference pattern along the waveguide axis. 

4 Conclusions 

We have considered the propagation of light in a planar waveguide ruled by the two-dimensional 
Helmholtz equation. Two Helmholtz equations with different refractive indices (3) form a super- 
symmetric pair in the sense of SUSY QM. 

Supersymmetry also describes the propagation of light beams of different colors in the same 
waveguide when the transversal index profile is elliptic. This profile determines the difference 
of light frequencies. SUSY is exact in the paraxial approximation, and broken beyond. In the 
paraxial regime, supersymmetric light beams of different colors have the same wavelength along 
the axis of the waveguide, giving rise to a stable interference pattern. 
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